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Holes from Faint X-ray Counts 
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ABSTRACT 

We investigate how hierarchical models for the co-evolution of the massive black hole 
(MBH) and AGN population can reproduce the observed faint X-ray counts. We find 
that the main variable influencing the theoretical predictions is the Eddington ratio 
of accreting sources. We compare three different models proposed for the evolution 
of AGN Eddington ratio, f^dd- constant /Edd — 1, /Edd decreasing with redshift, 
and /Edd depending on the AGN luminosity, as suggested by simulations of galactic 
mergers including MBHs and AGN feedback. We follow the full assembly of MBHs 
and host halos from early times to the present in a ACDM cosmology. AGN activity is 
triggered by halo major mergers and MBHs accrete mass until they satisfy the observed 
correlation with velocity dispersion. We find that all three models can reproduce fairly 
well the total faint X-ray counts. The redshift distribution is however poorly matched 
in the first two models. The Eddington ratios suggested by merger simulations predicts 
no turn-off of the faint end of the AGN optical luminosity function at redshifts z > 1, 
down to very low luminosity. 

Key words: cosmology: theory - black holes - galaxies: evolution - quasars: general 



1 INTRODUCTION 

Several hierarchical models for 
MBH and AGN populat ions 
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Volonteri et al 

2006) have proved successful in reproducing the AGN 
optical luminosity function (OLF) in a large redshift range 
(1 i; z ^ 6). Typically, these models assume that AGN 
activity is triggered by major mergers. Galactic interactions 
trigger gas inflows, and the cold gas may be eventually 
driven into the very inner regions, fueling an accretion 
episode and the growth of the nuclear MBH. Hydrodynamic 
simulations of major mergers have shown that a significant 
fraction of the gas in interact ing galaxies falls to the 
cente r of the merged system ( Mih os fc Hernauisd 11994 
1996): the cold gas may be eventually driven into the 
very inner regions, fueling an accretion episode and the 
growth of the nuclear BH. This last year has been especially 
exciting, as the first high resolution simulations of galactic 
mergers including BHs and AGN feedback ; showed that 
the merger scenario i s gene rally correct dSpringel et alJ 
2005 



nerge r scenario i s gene rally correct 
iDi Matteo et all 121)1)5 ). In hierarchical models of 
galaxy formation major mergers are responsible for forming 



bulges and elliptical galaxies. Support for merger driven 
activity therefore comes from the observed correlation 
between bulge luminosity - or stellar velocity dispersion - 
and black hole mass, suggesting a single mechanism for as- 
sembling black holes and forming spheroids in galaxy halos 
jMaeorrian et all 1993: iMarconi et all20oSlGebhardt et alJ 
l200CtlFerrarese fc Merritdl2000Tl . 

Notwithstanding the success at high redshift, hierarchi- 
cal models struggle to match the AGN OLF at low redshift 
(z < 1) by overpredicting the bright end, and underpre- 
dicting the faint end of the OLF, as the decrease of the halo 
merger rate with time is less dramatic than the observed fall 
of the AGN population. The overprediction of bright AGN 
can be imputed to inefficient cooling in large halos. Impos- 
ing an upper limit to t h e AGN host halo mas s (~ 10 13 ' 5 Mq, 
IWvithe fc Loebl (I2003D : iMarulli et alJ J2006M in fact signif- 
icantly improves the match at the bright end of the OLF. 

The under-abundance of faint AGN can be instead at- 
tributed to the assumption, common to most models, of 
very efficient accretio n, at rates clo s e to t he Eddington rate. 
iMerloni et"ai] (l2003fi and iMerlonil (12004 have shown that 
low-redshift AGN are probably accreting inefficiently, i.e. 
both at an accretion rate much smaller than the Eddington 
rate and with a low radiative efficiency. These considerations 
suggest that successful predictions for the evolution of AGN 
luminosity should include more sophisticated models for ac- 
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cretion. A first step in this direction can be taken by con- 
sidering the results of the recent merger simulations, which 
track also accretion on a central MBH. Although these simu- 
lations lack the necessary resolution for resolving the accre- 
tion process in the v icinity of the MBH, empirical models 
(Hop kins et al-lfeoOSl) . based on coupling results from the 
above simulations with the observed LF in the hard X-ray 
band (HXLF), have been shown to reproduce simulaneously 
severa l optical and X-ray observations. The iHopkins et alJ 
( 2005) empirical models, however, are not embedded in a 
cosmological evolutionary framework, that is they derive 
the MBH population properties at a given time, but not 
how the population of black holes at an earlier time evolves 
into t he MBHs present at a later time. From the observed 
HXLF IHopkins et alJ J2005I) derive the rate at which AGN 
of a given luminosity at the peak of activity must be cre- 
ated. This information is then used as a proxy for the galaxy 
merger rate which should provide the boundary conditions 
for determining the evolution of the MBH p opulation. 

W e here couple the predictions from Hop kins et alJ 
(2005) with the merger rate expected in the currently 
favoured cold dark matter scenario. The main novelties 
of the present investigation with respect to IHopkins et alJ 
(2005), are therefore that (i) the rate of mergers, which 
trigger AGN activity, is directly derived in the currently 
favoured cold dark matter (CDM) co smology. In pr i nciple , 
the empirical merger rate derived by IHopkins et alJ (2005) 
is not granted to correspond to the CDM one. Also, (ii), we 
grow MBHs in a self-consistent way, that is we trace the 
whole accretion history of MBHs from early times to the 
present, requiring continuity in the population. 

Our aim is to investigate here to which extent hierar- 
chical models, coupled with the prescriptions based on the 
above simulations (§1), can reproduce the low-redshift evo- 
lution of AGN. We show that reproducing the HXLF is a 
necessary, but not sufficient, condition for matching the red- 
shift distribution of faint X-ray counts (§2). We identify the 
redshift distribution of faint X-ray counts as the most sen- 
sitive observational result to discriminate between models 
(§3, §4). Finally, in §5 we summarize the results and discuss 
their implications. 



2 FORMATION OF MASSIVE BLACK HOLES 
AND GROWTH BY MASS ACCRETION 

In our framework pregalactic 'seed' holes form at early 
times. In mo st of our calculations we follow fVolonteri et alJ 
( 2003. l2005h . assuming that seed MBHs form with inter- 
mediate masses (m scc d ^ 600 Mq) in halos collapsing at 
z = 20 from rare 3.5-a peaks of the primordial density field 
dMadau fc Reesll200lt) as end-product of the very first gen- 
eration of stars. The assumed 'bias' assures that almost all 
halos above 10 11 Mq actually host a BH at all epochs. We 
also check the influence of the i nitial conditions by consid- 
ering seed MBH formation as in iKoushiappas et alJ ll2004l) . 
In this model seed MBH form from the low angular mo- 
mentum tail of material in halos with efficient gas cool- 
ing. In first approximation, seed MBH form in halos with 
mass above the threshold M H ^ 10 7 M Q (1 + z/18)" 3/2 , 
with a mass m Beed ^ 5 x 1 4 M Q (M H /10 7 M )(l + z/18) 3/2 
dKoushiappas et al.ll2004T) . We have dropped here the de- 



pendency on the halo spin parameter and gas fraction, as 
we are not interested in the detailed seed formation process, 
but only in testing an alternative model for seed formation 
which predicts much larger seed masses. 

Nuclear activity is triggered by halo mergers: in each 
major merger the hole in the more massive halo accretes 
gas until its mass scales with the fifth power of the circular 
velocity of the host halo with a normalization which repro- 
duces the observed local correlation between MBH mass and 
velocity dispersion (mBH — c* relation). 

The rate at which mass is accreted scales with the Ed- 
dington rate for the MBH. In model I the accretion rate is 
set exactly to be the Eddington rate. Defining the Eddington 
ratio as /Edd = MJMsdd, model I has a constant /Edd = 1. 
IShankar et alJ (l2004f) suggest that if the Eddington ratio 
evolves with redshift the MBH mass function derived from 
a deconvolution of t he AGN LF agrees better with the local 
MBH mass function ( Aller & Richstonc 2002: ]Shankar et alJ 
12004 iMarconi et all2004MShankar et al.N2004h suggest the 
following parameterization: 

f , \ J /Edd,0 z > 3 

/Edd(*) = j fEddA{1+z)/4] ^ z<3 (1) 

with /Edd,o = 0.3. Face value, Equation Qunderpredicts the 
LFs of AGN at high redshift, in our framework. As we start 
from small high redshift seeds, for our model II we modify 
EquationQas follows: 

r 1 2^6 

/ E dd(z) = < 0.078(1 + z) 2 - 0.623(1 + z) + 1.545 3 S z < 6 
I /Edd,o[(l + 2)/4] 14 z<3, 

where the quadratic form smoothly joins the z < 3 and z > 6 
functional forms. 

Finally, we test models where the Eddington ratio scales 
with the AGN luminosi ty. In model Ilia we p arameterize the 
accretion rate following lHopkin s et alJ ||2005| ). The results of 
simulations are presented in lHopkmsetH!uT<l2005li in terms 
of AGN luminosity. As our main variable is the black hole 
mass, we introduce some simplific ations to the model, but 
the general trend is preserved (see lHopkins et ail J2006t) for 
a thorough discussion). The time spent by a given AGN 
per lo garithmic interval is approximated bv lHopkins et al] 
ll2005h as: 

^ = |ateL_1 (lofe) ' (3) 

where t Q ~ 10 9 yr, and a = -0.95 + 0.32 log(L pcak /10 12 L Q ). 
Here I/ pe ak is the luminosity of the AGN at the peak of 
its activity. IHopkins et alJ (|200rt ) show that approximating 
ipoak with the Eddington luminosity of the MBH at its fi- 
nal mass (i.e., when it sets on the msH — c* relation) the 
difference in their results is very small. If we write the ac- 
cretion rate in terms of the time-varying Eddington rate: 
M = /Edd(£)M Ed d = /Edd(i)m B H/tEdd (*Edd = 0.45 Gyr), 
the AGN luminosity can be written as L — e/Edd(i)AfeddC 2 , 
where e is the radiative efficiency 1 . Differentiating with re- 

1 We determine the radiative efficiency self-consistently track- 
ing the evolution of b lack hole spins throughout our calculations 
(Volonteri ct al. 2005). We adopt an upper limit to the radiative 
efficiency of e = 0.16, as this corresponds, adopting the standard 
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spect to the Eddington ratio, we can write a simple differ- 
ential equation for /Bdd(t): 

ffigddft) = /Edd ( t ) ( eM E ddC 2 \ a 

dt \a\t Q \ 1O 9 L J ■ 1 ' 

Solving this equation gives us the instantaneous Edding- 
ton ratio for a given MBH at a given time, and we self- 
consistently grow the MBH mass: 

M(t + At) = M(t) exp ( f JLf Edd ( t )—\ . (5) 
\J At tEdd e J 

We define a lower limit to the Eddington ratio /Edd = 10 -3 . 

As will be discussed in Section §4, we also consider 
a modification of model Ilia, where we include a much 
stronger dependence of /Edd on the galaxy velocity disper- 
sion, in practice we modify the exponent a in Equation [3] as 
follows: 




Although not physically motivated, as Equation |3 is not 
either, Equation El was inspired by th e trend in accretion 
rates shown bv lDi Matteo et all l|2005ft . The Eddington ra- 
tios found with this modification (model Illb) are typically 
lower than in model Ilia. Again we set /Edd > 10~ 3 . Model 
Illb is therefore representative of a simple attempt to de- 
crease further the typical accretion rate of MBHs at low 
redshift. 

The main assumptions regarding the dynamical evolu- 
tion of the MBH population in our models can be fo und in 
IVolonteri et"aT1 (I2003I. I2005F ) : IVolonteri fc Reesl ll200d) . 



3 LUMINOSITY FUNCTIONS 

We have calculated the luminosity functions by implement- 
ing the different accretion models within a comprehensive 
model for black holes evolution in a Cold Dark Matter uni- 
verse. The history of dark matter halos and their associated 
black holes is traced by merger trees (Volonteri, Haardt & 
Madau 2003). The evolution of the massive black hole pop- 
ulation traces the accretion and dynamical processes involv- 
ing black holes. We have assumed that accretion onto nuclear 
black holes is triggered by halo mergers, and we have then 
computed the accretion rate and luminosity of the active 
systems as described in the previous section. At every step 
of the simulations we apply the appropriate Eddington rate 
to accreting MBHs. For model I /Edd = 1 for all MBHs at 
all times. For model II, /Edd is only redshift dependent (see 
Eq.^, while for models Ilia and Illb we determine /Edd as 
a function of the black hole mass at the beginning of the 
timestep, and of the host velocity dispersion. 

The luminosity functions are computed selecting the 
black holes which are active at the chosen output redshifts 
[z — 0.5,1,2,3), and weighting each of them according 
to the Press & Schechter function. We derive the AGN 



bolometric luminosity as L = e/Edd(£)Afed dC 2 . We apply 
the bo lometric corrections and spectrum of iMarconi et all 
( 2004) to model the SED in the blue band. The spectrum 
of unabsorbed (here dubbed Type I) AGN is described 
by a power-law with photon index T = —1.9, exponen- 
tially cut-off at E c = 500 keV. The averaged SED of ab- 
sorbed (Type II) sources (i.e., sources with absorbing col- 
umn log (Nu/cm -2 ) > 22) is described by the same Type 
I spectrum for E > 30 keV, and by a power -law (continu- 
ously matched) with photon index V = —0.2 llSazonov et alJ 
2004) at lower energies. The Type II/Type I ratio is, in gen- 
eral, functio n of luminosity and reds hift. Here, we adopt the 
model #4 of lLa Franca et all i2005l) . which explicity allows 
for the redshift and luminosity evolution of the Nh distri- 
bution, providing the be st fit model to th e HXLF of the 
HELLAS2XMM sample llFiore et al.ll2003h . Error bars, at 
1-cr, for the theoretical LFs have been computed assuming 
Poisson statistics. 

We have compa red our theoretical O LF at different red- 
shifts to the OLF bv lCroom et alJ (I2004T) obtained by merg- 
ing the 2dF QSO Redshift Survey (2QZ), with the 6dF QSO 
Redshift Survey (6QZ). Figurcs[3aiid|2|show a n absorption- 
correc ted OLF against the best fit models bv lCroom et alJ 
i2004h . In order to guide the eye, we have extrapolated the 
OLF at the faint end bel ow the observational limit (dashed 
lines) adopting the same ICroom et all i2004h fit. We note, 
however that the theoretical OLF fl attens towards lower 
luminosities, as has been observed jRichards et al] 120051 : 
iHunt et al.ll2004l) . 

The agreement between the theoretical and observed 
OLFs is good for all models (I, II, III), which result almost 
undistinguishable in the observed luminosity range (solid 
line in Figures and |5J. Model I overpredicts the bright-end 
of the OLF at low redshift (z < 1) more substantially than 
models II and III; on the other hand models with subEdding- 
ton accretion tend to underpredict the bright end at high 
redshift. In fact, we cannot form MBHs massive enough to 
power bright AGN by z = 3 if models II or III are considered 
at all times 2 . 

The OLF is known to be a biased tracer of the accre- 
tion history of MBHs, missing the vast majority of Type II 
objectes. Moreover it sp ans a smaller ran ge in luminosity 
compared to the HXLF jUeda et al]|2003T) . 

When we compare the theoretical and observed HXLFs 
of unabsorbed AGN (Figures [3] and , large differences at 
the faint end become apparent. Model I largely underesti- 
mates the faint end at z < 1, while at z > 2 it agrees 
very well in the luminosity range probed by current sur- 
veys. Model II underestimates the normalization of HXLF 
at all redshift, although the shape is satisfactorily matched. 
In model II accretion is simply not enough to grow black 
holes massive enough to account for the bright end of the 
HXLF. If the normalization is changed from, / E dd,o = 0.3 to 
/Edd,o = 1 (cfr. Lapi et al. 2006), the model fares much bet- 
ter at high-z (basically corresponding to model I at z > 2), 
but incurs in the same issues of model I at lower-z. 



conversion for accretion from a thin disc, to a maximum spin 
parameter of the BH a = 0.9. This value was chosen in agree- 
ment with lGammie et al] 120041) simulations, which suggest that 
the maximum spin MBHs can achieve by coupling with discs in 
magneto-hydrodynamical simulations is a ~ 0.9. 



2 We have modified model III, assuming / E dd=l at 2 > 12 
in order to obtai n signi ficant growth of MBHs before z=6. See 
IVolont eri & Recs] 120061) for a discussion on the constraints set 
by the OLF of bright AGN at z = 6. 
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Figure 1. Luminosity function of AGN in the B-band corrected 
for absorption. Clockwise: z=0.5, 1, 2, 3. Green squares show 
model I (/Edd=l)i magenta pentagons model Ilia (/Edd luminos- 
ity dependent). Solid lines : 2QZ/6QZ LF. The dashed lines show 
the extrapolation to faint magnitudes. 




log(V L Bo) 

Figure 2. Luminosity function of AGN in the B-band corrected 
for absorption. Clockwise: z=0.5, 1, 2, 3. Blue triangles show 
model II (/Edd redshift dependent), red circles model Illb (/Edd 
luminosity dependent). 



Models Ilia and Illb fare better in reproducing the low 
redshift HXLF, but predicts a large population of faint AGN 
at z > 2 and slightly underestimate the bright end at high 
redshift (z > 2). We note here th at our approach differs 
from that bv IlL opjrins et alJ i2005l) . The starting point of 
iHopkins et alJ 120051) is the HXLF. from which they derive 
the quasar birth rate, and consequently the OLF and other 
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Figure 3. Luminosity function of AGN in the hard X-ray band [2- 
10 Kev]. Symbols as in Figure 1. The dashed lines show the Ueda 
et al. (2003) HXLF in the observationally constrained luminosity 
range. 

diagnostics. Our approach instead follows the evolutionary 
path of MBHs and AGN, that is the population evolves self- 
consistently along the cosmic epochs, according to the ac- 
cretion properties stated in §2. The HXLF becomes there- 
fore a constraint, ra ther than an input of the model as in 
IHopkins et ail(l2005l) . 

We have investigated the impact of the initial condi- 
tions, by applying model Ilia t o a scenario in which seed s 
are much more massive, as in iKoushiappas et"al] i2004l) . 
The resulting luminosity functions are negligibly different 
with respect to the corresponding models assuming smaller 
seeds. This is because observable properties are mainly de- 
termined by the accretion history rather than by the initial 
conditions. Differences arise only at luminosity around 10 42 
erg s _1 , as a bump at 10 42 erg s _1 , where the MBHs with 
a mass around that of the initial seeds are clustered, and a 
sharp decrease faint- ward of 10 42 erg s _1 . 



4 FAINT X-RAY COUNTS 

The luminosity functions (OLF and HXLF) are the most 
sophisticated analysis of the evolution of AGN as a function 
of luminosity and redshift. On the other hand the available 
surveys do not probe yet the extreme faint end where the- 
oretical models mostly differ, except at very low redshift. 
Figures[j|]and2]show that theoretical models predictions, at 
2 > 0.5 branch off at luminosities not yet sampled by the 
HXLF. Number counts are the results of integrating over in- 
trinsic luminosity and distance. Number counts are a weaker 
test than the luminosity function, as AGN with a wide range 
of intrinsic luminosities are included at each flux. Neverthe- 
less, they are the most direct probe of the AGN population. 
They are independent of cosmology, allow to probe further 
the faint population, where the HXLF is still prohibitive be- 
cause of spectroscopic flux limits. We therefore compute the 
expected X-ray counts for the same AGN population that 
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Figure 4. Luminosity function of AGN in the hard X-ray band [2- 
10 Kev]. Symbols as in Figure 2. The dashed lines show the Ueda 
et al. (2003) HXLF in the observationally constrained luminosity 
range. 

we used to determine the luminosity functions and compare 
the model results to the most recent determinations of X-ray 
counts and their redshift distribution. 




log(S) [erg s" 1 cm -8 ] 

Figure 5. Predicted log(iV)/ log(S') in the observed soft-X band 
[0.5—2 keV] for the different models. Dotted lines show the con- 
tribution of sources with z < 2, whereas dashed lines the con- 
tribution of sources with z > 2. Solid line is the sum of the two 
components. Points report data obtained with Chandra (dots; 
Moretti et al. 2003) and XMM (squares; Baldi et al. 2002), and 
the bow-tie indicate s the result of the fluctuation analysis of the 
Chandra deep field jBauer et al.ll2004) . 



4.1 Basic Equations 

The number of sources (per unit solid angle) seen in the flux 
range S, S + dS by an observer located at zq, is 



dN 
dttdS" 



(7) 



where dV c /dzdQ, is the comoving volume element per unit 
redshift per unit solid angle, and n c (z,S) is the comoving 
density of sources at redshift z, with observed flux in the 
range [S, S + dS] . The integrated flux of a source observed 
at zo is given by 



S : 



2 , L„(M)dv, 

4-Kdl{z ,z) ' 



(8) 



where v — i^o(l + z)/(l + zq), dr,(zo,z) is the luminosity 
distance between redshift zo and z, L V (M) is the specific 
luminosity averaged over the source lifetime (assumed to be 
only a function of the BH mass, M), and Ais is the rest-frame 
frequency bandwidth. 

The background specific intensity J„ (20) observed at 
redshift zo at frequency uq, is 



. , , (i+zo) 3 r , ,di 

J -oUo) = — / e v {z) — dz, 

J ZQ 



(9) 



where dl/dz is the proper line element, and the comoving 
specific emissivity e^(z) is given by 



dn c 



e„(t) = / dM I L„(t-t',M)-^dt 



dM tL v {M) 



dn c 
dtdM' 



(10) 



The second approximated equality holds once we consider 
the source light curve averaged over the typical source life- 
time r, assuming the formation rate of sources per unit mass 
as constant over such timescale. 



4.2 Number counts 

We computed the number counts predicted by our differ- 
ent assumptions concerning the evolution of the Eddington- 
parameter along the cosmic history. The soft and hard X-ray 
logN/logS are shown in fig. |5] and fig. EJ respectively. The 
total counts (solid lines) are divided into the contribution 
of sources at z < 2 (dotted lines), and z > 2 (dashed lines). 
Model results are compar ed to a compilation o f X-ray data 
from the Chandra (dots. iMoretti et alJ ll2003l) ). and XMM 
(sauares: lBaldi et alJ (120021) ) deep field surveys. The bow-tie 
indicates r esults of the fluctu ation analysis of the Chandra 
deep field feauer et al.ll2004) . 

Model I fails to reproduce the slope of the observed 
logN/logS, falling short in the number of bright objects, 
and slightly overpredicting faint AGN. Moreover, the counts 
are dominated by high redshift sources for fluxes below 
log S < -14.6. In Model II, the redshift distribution of AGN 
is somewhat pushed towards lower redshift, because of the 
relatively longer accretion time involved. The model under- 
predicts the counts both in the soft and hard X-ray bands, 
as BHs do not have enough time to grow. Model Ilia matches 
well the observed logN/logS in the soft band, but underpre- 
dicts the counts in the hard band. Finally, Model Illb gives a 
reasonable good description in both bands, though it slightly 
overpredicts counts at very faint fluxes. AGN number counts 
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log(S) [erg s" 1 cm -8 ] 

Figure 6. Predicted log(AT)/ log(S) in the observed hard-X band 
[2-10 keV] for the different models. Lines and points are the same 
of Fig. El 



are dominated, in the entire observed flux range, by low red- 
shift objects, the contribution of sources at z > 2 becoming 
significant only at fluxes as faint as the limits of the more 
recent surveys. 



4.3 Redshift distribution of X— ray selected AGN 

Aiming at constraining further the 4 different models em- 
ployed, we compare the predicted redshift distributions to 
the results of the Serendipitous Extragalactic X-ray Source 
Identification (SEXSI) program, a survey designed to re- 
solve a l arge fraction of th e 2-10 keV cosmic X-ray back- 
ground fcckart et alJlioOfil) . The survey covers 1 deg 2 for 



fluxes > 1 x 10 



14 ergs" 1 

erg s" 1 cm -2 . Given the large survey area, 



cm , and 2 deg for fluxes 



> 3 x 10" 

the SEXSI program minimizes the effects of cosmic variance. 
The catalog contains a total of 477 spect ra, among which 43 8 
have redshift and optical identification feckart et alJ l2006 ) . 
The Type I AGN redshift distribution of the SEXSI pro- 
gram and of our selected models are shown in fig. [7| Note 
that we splitted the original SEXSI data in order to match 
our definition of Type I/Type II sources. Moreover, we have 
convolved the predicted number counts wit h the sky cov- 
erage of the survey for different flux limits llHarrison et alJ 

Model I fails completely to reproduce the observed red- 
shift distribution. In particular, the model largely underes- 
timates the number of sources at z < 1.5: the distribution 
peaks at redshift higher than observed. A better match to 
data is achieved by Model II. The general shape of the dis- 
tribution is reproduced, although the model largely under- 
estimates the total number of sources, as already pointed 
out. Model Ilia and Illb are in reasonable agreement with 
the data. Model Ilia falls short to the data at z < 1, while 
Model Illb overpredicts the number of sources observed in 
the range 0.5 < z < 1.5. In conclusion, the best agreement 




3 1 
redshift 

Figure 7. Distribution of Type I AGN as function of redshift for 
different models (points) compared to the result of Eckart et al. 
(2005; solid line) from the SEXSI program. In order to compare 
our results to the observed distribution, we have convolved the 
predicted number c ounts with the sky co verage of the survey at 
different flux limits iHarrison et alj|2003t) . and we have split the 
data so to match our definition of Type I/Type II sources. 



with the observed logN/logS, and with the redshift distri- 
bution of source s is found assuming a luminosity dependent 
Eddington rate iHopkins et alJl2005f) . 

We also tested that, in a model in which BH seeds are 
more massive, as in lKoushiappas et alJ <l2004l) . the resulting 
log N / log 5* and redshift distribution do not differ signifi- 
cantly with respect to models with earlier, smaller seeds. In 
conclusions, different formation scenarios for seed BHs are 
difficult to discriminate on the basis of available X-ray deep 
field surveys. 



4.4 Unresolved X— ray Background 

According to lMoretti et alJ feOOSf) . the intensity of the to- 
tal X-ray background (XRB) is 7.53 ± 0.35 x 10~ 12 and 
2.02 ± 0.11 x 10" 11 erg s" 1 cm" 2 deg" 2 in the 0.5-2 keV, 
and 2-10 keV energy bands, respectively. A large fraction, ~ 
94%, of the soft XRB (SXRB) has been attributed to sources 
with fluxes exceeding 2.4x 10~ 17 erg s" 1 cm -2 , while ~ 89% 
of the hard XRB (HXRB) is resolved into sources whose 
flux is ^ 2.1 x 10~ 16 erg s" 1 cm" 2 ( Moret ti et al. 2003). 
More recently, Hick ox fc Markevitchl feOOfj) estimated the 
unaccounted fraction of the XRB due to extragalactic un- 
resolved sources as 1.77 ± 0.31 x 10" 12 erg s" 1 cm" 2 deg" 2 
in the soft X-ray energy band (0.5-2 keV) and 3.4 ± 1.7 x 



12 erg s 1 cm 2 



deg in the hard X-ray energy band 



10" 
(2-8 keV). 

Using our different models, we compute the contribu- 
tion to the unresolved XRB due to faint AGN lying below 
sensitivity limits of current X-ray surveys. The cumulative 
contribution from sources with flux above a given thresh- 
old is shown in Fig. |H| where it is compared to the recent 
estimate of Hickox & Markevitch (2005; shaded area). Dif- 
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Figure 8. Predicted cumulative contribution to the unresolved 
XRB from different models as function of the X— ray flux. Different 
lines refer to different model: Model I (long-dashed line), Model 
II (dotted line), Model Ilia (dot-dashed line), and Model Illb 
(short-dashed line). The shaded area show the measured unac- 
counted background as reported by Hickox & Markevitch (2005). 
Top panel: unaccounted XRB in the observed soft— X band [0.5-2 
keV]. Bottom panel: unaccounted XRB in the observed hard— X 
band [2-8 keV]. 



ferent line styles refer to different models: Model I (long- 
dashed line), Model II (dotted line), Model Ilia (dot-dashed 
line), and Model Illb (short-dashed line). The unresolved 
SXRB and HXRB are shown in the top and bottom panel, 
respectively. Although significant differences (within a fac- 
tor of ~ 2) are found between different models, all pre- 
dict a contribution to the unresolved XRB consistent with 
available limits. Note that all models can account for the 
whole unresolved HXRB, while they give at most 50% of 
the unresolved SXRB. Our results imply the existence of 
a further population of faint X-ray sources in the soft 
band. Indeed, Salvaterra et al. (2006) have found that a 
significant contribution to the unresolved XRB may come 
from accreting BHs at very high redshift (z > 6), a red- 
shift range not considered here. By means of a dedicated 
model of the SMBH assembly at early times, consistent 
with the SDSS OLF at z — 6 and with ultra-deep X-ray 
constraints, they found that the contribution to the SXRB 
of very high redshift, undetected AGN could be as high as 
~ 0.4 x 10~ 12 erg s _1 cm -2 deg -2 , providing the residual 
unresolved flux. The contribution of such population to the 
HXRB is still consistent with the available limits, being only 

deg" 2 . 



0.9 x 10" 12 erg s _1 cm" 2 



5 DISCUSSION 

We have attempted in this paper to place constraints on the 
global accretion properties of the MBH population at z < 3. 
We consider the full cosmological evolution of MBH embed- 
ded in their host halos, rather than adopting an empirical 



approach which takes as a starting point the observed LF in 
a given band in order to explain the properties of AGN in 
other bands. We focus here on the strength of accretion, pa- 
rameterizing the accretion rate as a function of the Edding- 
ton rate, /Edd- We show that simple models which assume 
/Edd = 1 (model I), although highly idealized, still are able 
to explain satisfactorily the growth of MBH, as traced by 
the OLF, in a large redshift range, as claimed previously by 
various investigations. Even at low redshift, where hierarchi- 
cal models start to fail in their predictions, the OLF is very 
well reproduced. The HXLF probes fainter sources than the 
OLF, and the simplistic model I is shown to provide a poor 
match with t he HXLF at z < 1. A decreasing average Ed- 
dington ratio (IShankar et aljEo04l . model II) provides a bet- 
ter agreement with the shape of the HXLF, but underesti- 
mates the normalization, as there is not enough time to grow 
high mass MBHs if the Eddington ratio is not large at z > 3. 
A model ( model HI) with an E ddingto n ratio depend ing on 
luminosity (iHopkins et alJl2005l : Ishankar et alJl2004) seems 
to be the best match at low redshift. The difference among 
the models is magnified when the faint X-ray counts, and in 
particular their redshift distribution, are calculated. Model 

1 has a strictly hierarchical growth of MBHs, and moreover, 
a univocal relationship between MBH mass and AGN lumi- 
nosity. Model II allows massive black holes at low redshift 
to shine at lower luminosity, as the accretion rate, in units 
of the Eddington one, decreases with time. A rigid redshift 
dependence, however, creates two problems: first, subEd- 
dington accretion at z > 3 implies a much slower growth for 
high redshift AGN, thus causing an underestimate of the 
LFs normalization at all the considered redshifts. Second, 
the evolution of the Eddington rate is not fast enough, at 

2 < 3, to account for the faint end of the LFs, in particular 
the HLF. The accretion rates, and Eddington ratios, pre- 
dicted by simulations are here embedded into a cosmological 
framework (models Ilia and Illb). The resulting AGN pop- 
ulation provides a satisfactory match with the OLF, HXLF 
and faint X-ray counts at low redshift, while sources at high 
redshift are more problematic. The low accretion rates pre- 
dicted by simulations imply very long growth timescales for 
black holes, and therefore underestimate the occurrence of 
bright quasars powered by billion solar masses black holes 
at high redshift (z > 2). 

It is not clear, however, if the accretion rate found in 
simulations, which relates to model III, i.e., strongly subEd- 
dington for low-luminosity sources (and therefore in all cases 
for small MBHs) ind eed applies at very hi gh-redshift. The 
simulations on which IHopkins et alJ (120051) model is based 
assume mergers between "normal" galaxies. How "normal" 
are galaxies at z > 3? It might indeed be possible that the 
conditions in pre-galactic structures at high-redshift (e.g., 
the disturbed morphological state of galaxies) cannot be 
studied with simulations of mergers of evolved discs and 
bulges, but would require different initial conditions for the 
merging galaxies. Future deep surveys can help us distin- 
guish between various models in the high redshift Universe 
(see Figure|3J, and can probably locate the time, if any, for a 
transition between messy mergers with on average efficient 
accretion onto the MBHs, and standard galactic mergers, 
predicting long periods of inefficient accretion. 

Finally, we computed the contribution to the XRB of 
faint AGN lying below the sensitivity limits of current X- 



8 M. Volonteri, R. Salvaterra & F. Haardt 



ray surveys. We found that all model predict a contribution 
to the unresolved XRB consistent with the available limits, 
accounting for the whole unresolved XRB in the 2-8 keV 
band and for ~ 50% in the 0.5-2 keV. The residual back- 
ground intensity in the soft band may be provided by AGN 
shining at z > 6 (Salvaterra et al. 2006). 
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